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We present an experimentally realizable, simple mechanical system with linear interactions whose 
geometric nature leads to nontrivial, nonlinear dynamical equations. The equations of motion are 
derived and their ground state structures are analyzed. Selective "static" features of the model are 
examined in the context of nonlinear waves including rotobreathers and kink-like solitary waves. 
We also explore "dynamic" features of the model concerning the resonant transfer of energy and the 
role of moving intrinsic localized modes in the process. 
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I. INTRODUCTION. 

In the past few years, there has been a dramatic in- 
crease of interest in the behavior of soHtary waves and 
intrinsic locaUzed modes (ILMs) in dynamical lattices; 
see e.g., IJ for a number of recent reviews. One of the 
key reasons for this focus of interest has been the ability 
of such modes (which are ubiquitous in nonlinear lattice 
models) to localize the energy and transfer it in a tar- 
geted way 0, 0] ■ An additional, related feature of these 
modes is the important role they play in the conduction 
of heat (or equivalently transport of energy) along such 
dynamical simple lattices and how this relates to fun- 
damental macroscopic laws of thermodynamics such as 
Fourier's law of heat conduction; see e.g., Q for a recent 
review. Another application may be found in modelling 
the deformation and fracture behaviour of continuous 
structured media with internal degrees of freedom. Some 
initial model equations in this direction can be found in 

i 

On the other hand, an increasingly important theme 
in nonlinear physics concerns the interplay between non- 
linear dynamics and geometry, especially in lattice set- 
tings. The relevant contexts vary from long-range in- 
teractions on a fixed curved substrate to substrate- 
feedback models and coupled atomic chains 0, and 
from junctions between lattices with different masses 
to semi-circular polymer-like chains lld| and models of 



geometrically nontrivial DNA pd| . The unifying princi- 
ple in all these situations is that the geometry can sig- 
nificantly affect the static properties of excitations in the 
lattice (e.g., multi-stability), as well as dynamical ones 
(e.g., a variety of outcomes in the interaction of intrinsic 
localized modes with curvature). 



Motivated by these two emerging themes of nonlin- 
ear lattice dynamical systems, we propose in this work 
a mechanical example, in which even though the under- 
lying interactions of the system consist of linear springs, 
the geometry renders the interactions nonlinear. This, in 
turn, leads to complex features of the straightforwardly 
realizable (in a mechanical experiment) system at hand. 
Such examples are manifested in the coherent structures 
of the system rangin g fro m the familiar kinks to the more 
exotic rotobreathers [l^l (or other structures such as do- 
main walls or kinks with embedded defects |l3j|). They 
are also evident in the resonant energy transport features 
of the model that relate the energy conduction with the 
moving intrinsic localized modes (ILMs) that are present 
in the model (see below). We touch upon each of the 
above aspects to give a flavor of the rich and diverse 
properties that such a "geometrically-induced nonlinear- 
ity" model is endowed with. 
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II. 



MODEL 



The mechanical example that we examine consists of 
masses ("beads") that slide on fixed rings of radius R 
(even though the radius can be variable from ring to ring, 
we will here consider it to be fixed). Furthermore, the 
centers of the rings are at distance L between them. The 
chain of beads moving azimuthally along their respective 
rings is coupled through linear, elastic strings of a nat- 
ural length Iq; see Fig. ^^a) for a schematic. There are 
two interesting subcases. The rings can be in the same 
plane or they can be in different planes. While the latter 
case is also of interest, we restrict ourselves to the former 
(coplanar rings) in what follows. Then, 

Xn = nL + R cos 9,1, ?;„ = i?sin(6'„), (1) 

and the distance between adjacent particles is given by: 

(2) 



(yrn 



Vn) 



The model Hamiltonian for linear, elastic interactions 
reads: 



H 



E 



1 



-MR'' 



dt 



1 



(3) 



Out of the 5 model parameters (particle mass M , spring 
constant disk radius i?, distance between disk centers 
L and string natural length Iq), two (M, K) can be scaled 
out while out of the remaining three length scales, one 
can be used as a measure for the others (hence we set R — 
1 in what follows). The resulting equations of motion 
read: 



[sm(«„+i 
IqL sin 



'„) - sin(6l„ - 6n-i 
1 1 



sm(y„+i - ^n) sin(6'„ - 6'„_i)\ 
to 1 : ) , (4) 



^n — l.n 



with r2 „_^i = i2 + 2 - 2 cos(6l„+i - 6l„) + 2L[cos(6'„+i) - 
cos{6n)]- For convenience, we introduce the notation 

Ground States. The model of Eq. supports nu- 
merous complex structures (a number of which we will 
examine below), but there are only three types of the 
ground state structures bearing a very simple form. The 
parameter space {L,Iq) is naturally divided into three re- 
gions, depending on the corresponding types of ground 
state structure. We now summarize these structures (see 
also Fig. mb)). 

For Iq < L, we are in regime I, where the lowest energy 
structure is 

0n = c^, (5) 

for constant (f). Note that from static considerations, 
structures with any <j) are in indifferent equilibrium, but 



in the presence of dynamic perturbations, only structures 
with = 0, TT are stable. The physical origin of the insta- 
bility for (j) ^ {0, tt} can be easily understood considering 
the situation with all particles being at rest with 9n = (f) 
and one particle oscillating with a finite amplitude near 
the equilibrium position. Taking into account higher or- 
der anharmonic terms, one can demonstrate an asymme- 
try in the torques acting on the particle with respect to 
deviation from the right versus deviation from the left. 
This asymmetry gives rise to the force driving the par- 
ticle toward the closest stable position, namely (j> — or 
(j) = TT, where the symmetry is restored. There are two 
more symmetric structures corresponding to (f> = 7r/2 and 
(f) — — 7r/2, but they are unstable. To illustrate this issue, 
for the chain of iV = 400 particles, we set the initial con- 
ditions 0„ = 0, 0„ = -|- r„, with different magnitudes 
of (p and r„ being a random number homogeneously dis- 
tributed on [—0.05,0.05]. In Fig. QJc) we plot the time 
evolution of < 6'„ >= J2n^n- One can see that 
there are two stable positions (0 = and = tt) with 
respect to dynamic fluctuations. 

We thus restrict our considerations to the stable equiv- 
alent structures, = 0, tt. In this case, the dispersion 
relation reads: 



2jl- ^sin(^A:). 



(6) 



In the limit Iq — L, the linear spectrum collapses to a 
single point (in this case w = 0). Such geometrically 
induced limits where the linear spectrum collapses to a 
single point for all k are an interesting feature of the 
present model (see below) and will be called purely an- 
harmonic (PA) limits, as the harmonic linear part of the 
spectrum is eliminated in this case. This is a situation 
that bears some resemblence to the anti-continuum limit 
of models of ILMs Q. 

Regime II is defined hy L < Iq < Lq, where structures 
of the form: 



±[(-l)"0 + ^7r] 



/2 



L2 



(7) 



have zero energy (all springs have their natural length); 
6 = 0,±1. In this case, the linear spectrum is of the 
form: 



co{k) = ^pQ~L^^L^ + (4-pQ)sm\7rk), (8) 

which for Iq ~ 2 becomes fc-independent (PA limit), 
Lu{k) = Ly/A — L^. The spectrum of Eq. vanishes 
at Iq — L which is the border between the type I and 
type II regimes. 

Finally, regime III consists of the natural lengths such 
that Lq < Iq, when the ground state structure is 



0„ = ±(-l)"-. 
In this case, the linear spectrum is given by 



(9) 



.,.).2^(l-|)w(I-,.)+|l-l. (10, 
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One can see that for Iq = Lq/A, the width of hnear spec- 
trum vanishes (PA Umit), uj{k) = L. 

Having discussed the ground states, we now give repre- 
sentative examples of the interesting nonUnear behavior 
that is possible in each of the regimes highlighted above. 



III. 



REGIME I 



Regime I supports an interesting rotobreather with one 
rotating particle. Such an example is shown in Fig. EJa), 
where the particle n = rotates with angular velocity 
w = 20 at L = 3 and Iq = I- We show the strobo- 
scopic picture of motion with intervals T/10 (where T is 
the period of the solution). Particles with positive and 
negative n are practically at rest at their respective po- 
sitions. The neighbors of n = particle can be at rest in 
the positions where the averaged (over a period) torque 
acting from the moving particle is equal to zero. For a 
large rotobreather frequency, the angular velocity of the 
n = particle is almost constant and, hence, the averag- 
ing over time can be substituted by averaging over angle 
Oq. Thus, an approximate rotobreather solution can be 
given as follows: 



On = Lot , for n = , 
= , for n > , 
±(7r - 9i) , for n < , 



(11) 



i.e., the 0th particle moves with constant angular velocity 
Lo while all other particles are at rest at the positions 
defined through the angle Oi , which is a root of 



27T 



(12) 



where Mi(0q, 9i) is the torque acting on the particle n = 
1 from the neighboring particles under the assumption 
that 02 — 01 and 6*0 is arbitrary. While the solutions of 
Eq. H12|) with 0i = 0, tt lead to unstable configurations, 
the root lying in the interval tt/2 < 0i < tt leads to a 
stable rotobreather which very accurately captures our 
numerical findings (see the caption of Fig. |2Ia)). 

Note that in the case of roto-breather presented in 
Fig. 121 a), dynamic fluctuations are absent. Introduc- 
tion of such fluctuations would result in the appearance 
of a weak forces driving the particles with n < to the 
position 0n = and particles with n > to the positions 
with 0n ~ TT [see Fig. ^c)]. An equilibrium state in this 
case would look similar to what is presented in Fig. [^a) 
only in the vicinity of the roto-breather, but more distant 
particles would approach the positions which are stable 
against dynamic fluctuations. 



and provide a numerical example of a highly localized, 
moving ILM. 

Kink solution. Using 



(-ir(| 



(13) 



where e„ <C tt and assuming that e„ varies slowly with 
n, we obtain from Eq. [up to cubic terms] 



2 (^n-l 



2e„ + En+l) + CSr, 



B 
L 

^ [(e„+i -f- Enf + (en + e„-i)^] , (14) 



with B 



C 



41 



and D 



8 /i lo{L' + l) 

3 ir~ 



duces to the </? -equation. 



In the continuum limit, Eq. 114() re- 



£tt = Be^^ + Ce- De^ 



(15) 



When C > and D > 0, the background potential has a 
double-well structure and Eq. (I15|l supports topological 
solitons (kinks and antikinks) (IJi] . When the kink width 
is much greater than L, the solution of Eq. ifT^ can be 
used to approximate the kink solution of the discrete Eq. 



(-1)" 




tanh [Q{nL — vt)] 



(16) 



where v < ^/B is the kink velocity and Q = 
^C/[2(w2 — B)]. We have verified that even for rela- 
tively small kink widths (i.e., moderate discreteness), Eq. 
()16|l approximates well the numerically obtained kink so- 
lutions [see Fig. I3b)]. 

We note in passing that in regime II there exist vari- 
ous different types of domains, as illustrated in Eq. j?)). 
Hence, there are possibilities to create additional kinks, 
connecting different steady states than the ones presented 
above. 

Moving ILM. An interesting example of a moving ILM 
is presented in Fig. Hfor Iq = 2 (PA limit), L = 1.2. The 
localized mode emerged from the local perturbation in- 
troduced to n = particle of the structure On = (— 1)"</' 
by setting 6*0 = + 7r/2 with zero initial velocities for all 
particles at t = 0. The asymmetry in the displacements 
of particles renders this structure less tractable analyti- 
cally. We will give a very simple analytical approximation 
for the ILMs (standing and moving) in regime HI where 
the displacements are symmetric. However, we believe 
that the mechanism of ILM propagation in the regime 
II is the same as in the regime HI, i.e., the relay-like 
resonant energy exchange between particles. 



IV. 



REGIME II 



REGIME III 



As two prototypical examples of the structures that 
can be supported in regime II, we derive a kink solution 



Standing and moving ILMs. Regime HI supports ILMs 
which, depending on model parameters, can be moving 



4 



or standing. We have also found that the existence of 
standing ILMs precludes the existence of moving ones 
and vice versa. To better understand this phenomenon 
we first solve an auxiliary problem of resonant energy 
transfer and then give some simple analytical estimates 
for standing and moving ILMs. The robustness of moving 
ILMs against their collisions is also verified numerically. 

Here we consider the PA limit, Iq — L^/A. If a particle 
is forced according to: 



(17) 



with ^ <JC TT, we are interested in the motion of its near- 
est neighbors, 0-i{t) = 9i{t) = ~ e{t) with £ <C tt, 
assuming that all other particles are at rest at their equi- 
librium positions. Retaining up to linear terms in e and 
cubic in A, we obtain 



uJo ^ L, (3 = 



e + ujIs = 13 si-n^ [ujt) + ^A^ sin^(tji) 
L{2-L^) 2-L2 



1 



with the particular solution 



e{t) 



I3A^ 



PA^ 



2(a;2 - 4lo^) 



cos{2ojt) 



4(c.2-c.2) 
7^3 



2^4 1 , , 
(18) 



sin(a;i) 



4(t^2 _ 9^2) 



sin(3wt), (19) 



which gives the first three resonance harmonics. When 
the forced particle oscillates with a nearly resonant fre- 
quency, the amplitude of its neighbors can grow signifi- 
cantly. 

Equation (|18|l does not take anharmonicity into ac- 
count. The latter effect was studied numerically, where 
we found that the anharmonicity is hard, i.e., the oscilla- 
tion frequency of the n = 1 particle grows as a function 
of amplitude, for L < L* k 1.65, and the situation is 
reversed for L > L*. 

The dynamics of the chain with one forced particle 
differs qualitatively for hard and soft anharmonicity. For 
L < L*, when the free particle has maximum ampli- 
tude, it oscillates in phase with the forced particle, and 
its amplitude excceeds A in the resonance regime while 
for L > L* it does not because the particles oscillate out 
of phase. In other words, efficient inter-particle energy 
exchange occurs only for hard anharmonicities. For the 
chain with soft anharmonicity, any local perturbation re- 
mains local. 

We carry out the following numerical experiments. We 
excite a single site according to Eq. (|17|l for times < 
t < T with T = 5000 and calculate the power of the 
energy source, W = E/t, where E is the total energy of 
the chain at t = r. The chain is long enough so that at 
t = r the perturbation produced by the forced particle 
has not reached the boundaries. 

We have found that W can be nonzero only for L < L*, 
regardless of the magnitudes of A and ui in Eq. (|17|l . The 



mechanism of the energy transfer is the emission of mov- 
ing ILMs. In Fig. Efa), we present W as functions of uj 
for L = 1 and for A — 0.2, 0.3, 0.4. It can be seen that 
the smaller A is, the narrower the window of W > 0. 
The figure also shows the distribution of energy in the 
chain &t t — t for L — 1, Iq — -/jo/4, and driving pa- 
rameters A = 0.4 and (b) uj = 0.992, (c) uj = 1.020, 
(d) Lo = 1.033. In (b) and (d) (edges of the window 
with Vl^ > 0), the forced particle emits ILMs periodically 
while in (c) (central part of the window) chaotically. For 
L > L*, there is no efficient energy exchange between 
particles and moving ILMs are not possible. Instead, sta- 
ble standing ILMs arise that are localised at the excited 
particle. An approximate solution can be expressed by 
a conventional perturbation method assuming that only 
one particle moves: 



9o(<) = - +Ai sin(wi) + A3 sm{3ujt), 

+ -B^?, A3 = „ -TT, 

° 4 ^' ^ -yiijjl + TJBAl' 



(20) 



where 



2 _ 4;o(L''+2) 
^0 - 



2. B 



lo 



3L, 



The moving ILM is (practically) localized at three par- 
ticles: 



aAi sin ( utt ± 



2tt 



: Ai sin {ojt) + A3 sin (30-1^) , 
9n+i = aAi sm lujt^ — 



(21) 



where upper and lower signs correspond to ILM moving 
in positive and negative directions, respectively. The am- 
plitude Ai of the ILM is a free parameter. The relation 
between Ai, the ILM frequency uj, and the third har- 
monic amplitude, ^3, is given by Eq. H20(l . Empirically 
we have found that a = 0.4 gives a good result over a 
wide range of ILM amplitudes (see Fig. |SJ). 

In Fig. ini we show an in-phase collision of two ILM 
defined by Eq. ^ with A = 0.7, a = 0.4. Model 
parameters are L = 1, Iq = L^/A (PA limit). Shown are 
the snapshots of Hn at different time t, where Hn is the 
total energy of nth particle (kinetic and potential) . We 
can see that such moving ILMs can interact with each 
other in a quasi-elastic fashion. 

Rotobreathers. Regime III also supports a rotobreather 
with one rotating particle. Here again, as in Sec. IIIII 
we assume that for a rotobreather with sufficiently large 
frequency the torque acting from the rotating particle 
on its nearest neighbors can be estimated by averaging 
over angle. In Fig. \7\ we show the (averaged over an- 
gle) torque acting from rotating (n = 0) particle on its 
nearest neighbor {n = 1), which is assumed to be at rest 
at 9i, for lo = Lfj/4 (PA limit) and L = 2 (dotted), 
L — 1.5 (dashed), L = 1 (solid). We are interested in the 
positions 9i corresponding to zero torque. For L < 1.3 
there are only two roots (one of them is stable) and for 



5 



L > 1.3 there are four roots (two of them are stable). 
For L = 5, for example, the stable root was found at 
di = — 7r/2 + 0.084 which is in a good agreement with 
what is observed for the rotobreather presented in Fig. 
|S| For L — 5, another stable root for the n = 1 particle 
was found at ^±1 = 7r/2 + 0.135 and the existance of this 
rotobreather was also confirmed numerically (see Fig. |3- 

VI. CONCLUSIONS 

In this paper we have presented a novel nonlinear dy- 
namical system, consisting of an easily realizable mechan- 
ical example where the nonlinearity is induced by the ge- 
ometry of the problem. We have illustrated the laws of 



motion and the rich static, dynamic (both equilibrium 
and non-equilibrium) behavior of the system. We have 
identified some of the relevant coherent structures includ- 
ing kink-like heteroclinic connections and roto-breathing 
periodic orbits and have seen some of the interesting dy- 
namical phenomenology including the "conducting" (for 
hard anharmonicities) or "insulating" (for soft anhar- 
monicities) behavior of the system and the role of mov- 
ing or standing ILMs, respectively, as energy carriers. It 
would be of interest to examine further from an analytical 
(as well as from a numerical or experimental) perspective 
the phenomenological wealth of such a model. Such stud- 
ies are currently in progress and will be reported in future 
publications. 
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FIG. 1: (a) The bead-ring configuration: the rings of radii 
R arc at distance L between them (either on the same or on 
different planes). Tfic dynamical variable of interest for each 
particle is its azimuthal angle 9n and the chain has a nearest 
neighbor coupling through linear, elastic springs of natural 
length Zq. The distance between adjacent particles is denoted 
by rn,n+i. (b) The parameter space of the model {L,lo) di- 
vided into three regions with different ground state structures, 
(c) Illustration of dynamical instability in the regime I for all 
(p except for (p = and (p = ±7r. Average atomic positions 
< On > for the chain N = 400 particles are shown as the 
functions of time. Initially the particles are placed aX, On = <p 
with different magnitudes of (p, and a small amplitude random 
perturbation is introduced in the particle positions to initiate 
their vibrations. 
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FIG. 2: (a) Numerical solution for a rotobreather in regime 
I. The n = particle rotates with angular velocity oj = 20 at 
L = 3 and /o = 1- We show the stroboscopic picture of motion 
with intervals T/IO (a tenth of the period T). Particles with 
positive and negative n are practically at rest at 6l„ = 7r/2 ± 
0.261. Numerically we found the root of Eq. 1121 at Q\ — 
■k/2 + 0.259. (b) A kink solution is shown in regime II for 
V = 0.2, L = 1, = \/4:Sm^ + L'^, with (j) = (9/20)7r (close 
to <f> = n /2). Even and odd particles are shown by filled and 
open circles, respectively. Even in this case of (not very big) 
kink width, Eq. (1161 provides a very good approximation. 
The arrow shows the direction of propagation of the kink. 
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FIG. 3: Moving ILM in the regime II for Zo = 2 (PA limit) 
and L = 1.2. 



8 




FIG. 4: (a) The power of the energy source W is shown as 
a function of the driving frequency uj for L = 1 and for A = 
0.2, 0.3, 0.4. Right Paneb the distribution of energy ait — t. 
Model parameters are L = 1, lo = Lq/A, A = 0.4 and (b) 
UJ = 0.992, (c) UJ = 1.020, (d) uj = 1.033. The panel shows 
the particle energies E„, averaged over the period 2-k/u!. In 
(b) and (d) the forced particle emits ILMs periodically, while 
in (c) it emits chaotically. 
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FIG. 5: A moving ILM is shown by the functions 9n{t) for 
the four nearest nodes. Particles show the relay-hke motion 
oscillating near the equilibrium positions 6„ = ( — l)"(7r/2). 
Model parameters are lo = I/o/4 (PA limit in regime III), 
and 1/ = 1 (< L*). The parameters in the solution Eq. H21|l 
are A\ = 0.6, a — 0.4. Marginal radiation can be seen after 
the ILM passes a node (at large times). The ILM propagates 
rather slowly, it travels one lattice spacing L in about ST. 




FIG. 6: Elasticity of the in-phase collision of two moving 
ILMs defined by Eq. I^TJ with A = 0.7, a = 0.4. The model 
parameters are L = 1 (< L*) and lo ~ L'q/A (PA limit). 
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FIG. 7: Torque (averaged over the angle do) acting from ro- 
tating particle (n = 0) on its nearest neighbor (n = 1) which 
is at rest at Oi for lo = Lq/A (PA limit) and L = 2 (dotted), 
L = 1.5 (dashed), L = l (solid). 
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FIG. 8: Rotobreather in the regime III initiated by setting 
initial angular velocity 6o = 20 and initial positions of the 
n — ±1 particles 0±1 — —7r/2± 0.084 corresponding to zero 
averaged torque acting from n = particle. Model parame- 
ters arc L = 5 and lo — Lq/A (PA limit). Dashed horizontal 
lines show 9n = ±7r/2. 
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FIG. 9: Rotobreather corresponding to the root of equation 
< Ml >= which appeares at large L. In this case L — 5, 
lo = Lq/4: (pa limit) and the root was found for n — ±1 
particles at 7r/2 ±0.135. Note that here the particles n — ±1 
oscillate near n/2 but not near — 7r/2 as in the case presented 
in Fig. |H| This is because the stable root which appeares at 
L > 1.3 is shifted by, roughly, tt compared to the root existing 
for all L (see Fig. E}. 



